Mon. Not. R. Astron. Soc. OOP. [Tlflil C2012') Printed 12 April 2012 (MN I^TeX style file v2.2) 



(N- 

^ ■ 

o: 

a; 
< 



Evidence for inhomogeneous reionization in the local 
Universe from metal-poor globular cluster systems 

Lee R. Spitler^* Aaron J. Romanowsky^, Jiirg Diemand^, Jay Strader^'^, 
Duncan A. Forbes^, Ben Moore^ and Jean P. Brodie^ 

Centre for Astrophysics & Supercomputing, Swinburne University, Hawthorn, VIC 3122, Australia 
^ University of California Observatories, Santa Cruz, CA 95064, USA 
^ Institute for Theoretical Physics, University of Ziirich, CH-8057 Ziirich, Switzerland 
'^Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138, USA 
^ Menzel Fellow 

Accepted 2012 March 31. Received 2012 March 8; in original form 2011 September 8 



o 

o 

in: 

(N 

> 
m 

(N 
<N 



ABSTRACT 

Exploiting a fundamental characteristic of galaxy assembly in the ACDM paradigm, 
the observed spatial biasing and kinematics of metal-poor globular star clusters are 
used to constrain the local reionization epoch around individual galaxies. Selecting 
three galaxies located in different environments, the first attempt at constraining the 
environmental propagation of reionization in the local Universe is carried out. The 
joint constraint from the three galaxies {zreion = 10.5^Q g) agrees remarkably well 
with the latest WMAP constraint on Zreion for a simple instantaneous reionization 
model. More importantly, the range of Zreion values found here are consistent with 
the global range of zreion estimates from other observations. We furthermore find a 
1.7(7 indication that reionization completed in low-density environments before the 
intergalactic medium in high-density environments was reionized. This is consistent 
with certain theoretical models that predict that reionization was globally prolonged 
in duration, with neutral hydrogen pockets surviving in high-density environments, 
even after the surrounding regions were reionized. More generally, this work provides 
a useful constraint on the formation history of galaxy stellar halos. 

Key words: dark ages, reionization, first stars — galaxies: star clusters: gen- 
eral — galaxies: haloes — galaxies: kinematics and dynamics — galaxies: individual 
NGC 1407, Messier 87, Milky Way. 



1 INTRODUCTION 

Reionization marks the phase change of the Universe's in- 
tergalactic hydrogen gas from a neutral to an ionized state. 
Though the general reionization pictur e follows from the 
funda mentals of Big Bang cosmology (|Barkana fc Loebl 
l200lh . a detailed timeline of how it propagated throughout 
the Universe is difficult to constrain because most of the di- 
rect observational signatures of this event are > 12 billion 
light years distant from Earth. 

Despite this challenge, observations show that the in- 
tergalactic medium was reionized no later than a redshift of 

2 ~ 5, not before ^ ^ 12 and completed on a timesc ale of at 
least A2; > 0.06 jFan et al.ll2006l: lOuchi et al.ll2010l: Robert- 
son et al.T2010l: iBolton et al.ll2010l: iBowman fc Rogerdlioiol: 
iLarson et al.ll201lh . Though the dominant ionizing source 
contributing to reionization is still an open question, there 
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is general agreement that the ionizing sources first appeared 
in high-density environments (Diikstra ct al. 2004; F an et al.l 



20061: Lidzetal. 2007; Wise & Abel 2008; Power et al. 



Back et al. I2OIQ: Bunker et al.ii201Q: iRobertson et al 



Srbinovs kTfc Wvithell2010|: [Willott et aL[ [2010l: 
20101: Bouwcn s et al.ll20 11: Dopi ta et allboill: 



200' 



201 



Yan et al] 



Lorenzoni 



Srbi: 

201a ■ , 

et al. 2011; Mitra et al.l l2011. 201^ 

A key piece of information that is not yet constrained 
is the environmental propagation of reionization: did 
reionization complete in dense environments first or were 
low-density voids the first locations to ionize? While 
high-density environments of the Universe probably had 
a "head-start" and thus contained more ionizing sources, 
the same locations were also more likely to be shielded 
by dense clouds of dust and neutral hydrogen (HI) gas. 
Furthermore, if the sources were predominately small galax- 
ies, the UV background would have been more uniformly 
distributed than if the ionizing sources were typically active 
galaxy nuclei of rare, massive proto-galaxies. Also, unlike 



© 2012 RAS 



2 L. Spitler et al. 



ultraviolet radiation, X-rays produced by gas accretion 
onto black holes can escape local gas and dust absorption 
and "pre- heat" gas in lower-densitv regions fe.g. Baek 



et al. i2010f ). Indeed, large-scale radiative transfer models 
show that the ionizing front can travel vast distances, from 
massi ve galaxies in dense regions to low-density locales 
(e.g IWeinmann et al.ll2007l : iFinlator et aLll2009l : llhev et all 
I2OIII) . Understanding how reionization propagated through 
various environments will help constrain the properties 
of the ionizing sources and the state of the intergalactic 
medium during the reionization epoch. 

One possible surviving relic of the reionization epoch 
may be the bimodal distribution of globular star clusters 
withi r i large galaxies (Forbes et al.lll997t Ivan den Bergbl 



2001'; 'Beaslev et al."2002'; 



'Santos"2003'r 



Rhode ct al. 2005, 
2007, 2010; M oore et al. 2006; Bckki 2005; Bckki ct al. 2008; 
Griffen et al ll2O10l l: systems of glob ular cluster s show two 



distinct metalli citv subpo pulations (IZinnI Il985l: G ebhardt 
& Ki ssler-PatiglTg g 



19 99: Ncil sen fc Tsvetanovlll999l: Larsen 
et al. I2OOII: iK undu & Whitmore 2001; Peng ct al. 2006'; 

2006, 2007; Kundu & Zcpf 2007 ; S pitle r ct al. 



Strader et al.l ._ , , ^_ _ .^c:.^^^^, —- 

2008ah . The subpopulations include the "metal-poor" globu- 



lar clusters (MPGCs) with typical metallicities a few percent 
of Solar metallicity and the "metal-rich" globular clusters 
(MRGCs) with metallicities ~ 30% Solar. 

Various theoretical works have invoked reionization as 
a way to produce the bimodal metal l icity distribu t ion in 
GC syst ems (e.g. 'Beaslev et al."2002'; 'Santos^ '2003'; 'Bekki' 
2005; M oore et al . 2006; Bckki et al. 2008; Griffon ct al. 
201Q ) . More specifically, these studies explored models where 
MPGCs form within small (Myjr ~ lO^Mffi) dark matter ha- 
los at 2 > 6 (|Bromm fc Clarkdl20oj : iBolev et al]|2009l ). In 
this model, the MPGCs formed until cold molecular clouds 
within these halos was heated up during the reionization 
epoch. This suppressed globular cluster formation for a pe- 
riod of time and allowed gases within the host galaxy to 
enrich in metals. After some time had passed, a second gen- 
eration of globular clusters, the metal-ric h GCs, started to 



form during galaxy merger eve nts (e.g. lAshman fc Zep3 
I1992I : iBeaslev et al.l 12002: Griff en et al.l 120 10*) and/or the 
growth of galaxy disks ( i Shapiro et al.l boiOi) at redshifts 
2 ~ 2 - 4 (jShapiro et al.lboiGl : ISpitleij|201Gh . Although the 
global nature of the reionization epoch may provide a conve- 
nient way to explain the ubiquitous GC metallicity bimodal- 
ity observed in nearby galaxies, a number of interesting al- 
ternative models have been proposed to explai n GC metal- 
licity distributions without reionization (e.g. I Cote et al.l 
Il998 >: Scanna pieco et al.ll200 4': Pipino et al.ll2007l: Hasegawa 
et al. ,2009: .Shapiro et al.ll2010. : .Muratov fc GnedinI [2OI0I: 
ICrav fc Scannapiecol 20101 ). 



If reionization truncated globular cluster formation, the 
typical ages of metal-poor globular clusters should correlate 
with the local reionization epoch. This means metal-poor 
globular clusters can be used to measure ZrEion for indi- 
vidual galaxies. Current MPGC age estimates overlap with 
the expected reionization epoch. Galactic MPGCs are con- 
strained to be older than ~ 11 Gvrs or z ~ 2 ( Krauss fc 
Chabover BoOSi ). with more stringent constraints on individ- 
ual MPGCs suggestive o f a pre-reionization epoch formation 
(e.g. [Hansen et al.ll2004l. though exceptions do exist: Hansen 
et al. l2007f ). Constraints on the ages of individual extra- 



galactic MPGCs are poorer (> 10 Gvrs. see refs. in Brodie 
fc Strader l2006ir but are consistent with a pre-reionization 
formation. Also, within the 1-sigma uncertainties the mean 
MPGC metallicities of various early-type galaxies are statis- 
tically consistent with prel iminary measurements for galax- 
ies a t redshifts z ~ 6 — 8 JPinkelstein et al.ll2010l: Bouwens 
et al. I2OIOI : iLabbe et al.ll2010l ). This is shown in Fig. [1] and 
is consistent with a scenario where MPGCs were in place at 
redshifts greater than z ~ 7, well into the expected reion- 
ization epoch. 

Unfortunately, the limited observational information 
about extragalactic globular clusters and large modeling un- 
certainties mean that direct methods to age-date extragalac- 
tic globular clusters do not provide a useful constraint on the 
local reionization epoch around individual galaxies. 

A novel way to determine when the metal-poor GC for- 
mati on epoch finished was recently proposed bv Diemand 
et al. l|2005^ and lMoore et al.l l|2006l) . The technique is based 
upon a fundamental property of a ACDM universe, where 
the most massive or rarest dark matter halos (and anything 
that formed within them) will tend to be centrally concen- 
trated within their final host halos and show unique kine- 
matic signatures at the present epoch. Thus by constraining 
the spatial and kinematic properties of metal-poor globular 
clusters today, information about their progenitor halos can 
be recovered. Since the rarity of a halo for a given virial mass 
depends on redshift, this information can be used to age- 
date metal-poor GCs with enough accuracy to constrain the 
reionizatio n epoch around the galaxv hosting the GCs (Die- 

Bekk i 2005). 
Moore et al 



reionizatio n epocji around the galaxv hosting 
mand et al. l2005l : iMoore et al.ll2006l see also 

This technique was demonstrated in 

(I2OO6I) using the Milky Way's globular cluster system. They 
found evidence for a reionization red shift of Zre.ion = 10 ± 2 
(updated for a WMAP7 cosmology, iKomatsu et al.ll201lh . 
In the present work, this technique is applied to two 
additional galaxies located in denser environments than the 
Milky Way to help constrain the temporal propagation of 
reionization through the local Universe. 

Section [2] provides a detailed background, describing 
in detail both the theoretical framework of the "Diemand- 
Moore" methodology and relevant globular cluster observa- 
tions. In SectionO the observational data are presented. The 
main mass modeling and Zreion measurements are presented 
in Section (31 It is followed by a discussion in Section [S] 



2 METHOD AND ASSUMPTIONS 

It has been shown that the growth of structure in the Uni- 
verse can be thought of as a Gaussian mass-density field 
where m ass-density fluctuations form and assemble h ierar- 
chically (jPress fc Schechtedl 19741 : iLacev fc Colelll993D . The 
rare mass-density fluctuations in the Univers e are biased to- 
wards high-density regions o f the Universe (|Cole fc Kaiser! 
ll989l:ISheth fc Tormenlll999h . As a consequence, the rarest 
fluctuations will tend to be the most cen t rally-concentrated 
material in a galaxy l|Moore et al.ll 19981: IWhite fc Springell 
I2OOOI) . By convention, the height of the fluctuations are rep- 
resented by V, the number of standard deviations above the 
mean mass-density level at that epoch. 
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Figure 1. Mean mctallicity of metal-rich (circles) and metal- 
poor (squares) globular clusters in Virgo cluster galaxi es as a 
funct ion of their hos t galaxy stellar m ass (for details see ISpitled 
l20ld. GC data from lPeng etal]|2006l) . Metallicity measurements 
for galaxies at redshifts 2 ~ 6 — 8 are given with measurement 
uncer tainty ranges (blue triangle s; data from iFi nkelstei n et alj 
l201Ct see also the simulations of ISalvaterra et aL[l201 ll) . Using 
metallicity as a timescale indicator, the overlap in stellar metal- 
licity between metal-poor globular clusters (MPGCs) and the 
high-redshift galaxies supports t he idea that MPGCs for med at 
redshifts 2 ~ 7 — 10, when the iFinkelstein et al. I 1*2010) galax- 
ies formed their stars. Dashed lines are observed mass-metallicity 
(from io nized gas observations) relations for red shifts z ~ 2.2 and 
z ~ 3.5 llErb et ahlliooj: iMannucci et al.ll2009h . As discussed in 
ISpitlej \20idt). the over ap between these relations and the metal- 
rich glol )ular clusters suggests they formed later than the MPGCs, 
at 2 ~ ; — 4. 



Diemand et al. l|2005l) used N-body simulations to 
explicitly show that the location and kinematic properties 
of mass structures within a galaxy today depend on the 
rarity of the mass structure or "halo" when it reached 
its Jeans mass and gravitationally collapsed. They found 
that the rarest halos are more likely to accrete onto a 
galaxy at earlier times, when the gravitational potential 
of the galaxy was relatively small. Thus any object that 
formed within a rare halo is more likely to have a smaller 
infall velocity and will occupy the central regions of the 
galaxy until the present epoch. This means that an object's 
location and orbit or kinematics within a halo today can be 
used to recover the properties of the halo it formed within 
l|Diemand et all I2OO5I : iMoore et all I2OO6I ). Examples of 
"tracer" objects include a galaxy's satellites, its stellar halo, 
remnants of population III stars, and metal-poor globular 
clusters. 

The "Diemand-Moore" method is a way use the spa- 
tial concentration of MPGCs relative to the hosting galaxy's 



halo mass profile to constrain the rarity, va, of the halos that 
the MPGCs formed within. Results from computer simula- 
tions are used to translate the spatial bias of the MPGCs 
into a va constraint. The lya value alone does not constrain 
the redshift when MPGCs stopped forming, so an additional 
assumption is made to break the degeneracy between red- 
shift and the typical halo mass for a given ly. By invoking 
the argument that star formation only occurs efficiently in 
halos with virial temperature more than Tvir = lO" K, the 
degeneracy can be broken and the redshift when MPGCs 
were suppressed by reionization is found. More details are 
provided in the following Sections. 

The Diemand-Moore method requires at least three in- 
gredients to constrain the lya value of MPGC progenitor 
halos: MPGC surface density profile from observations, a 
halo mass model for the ga laxy ho sting the MPGCs and a 
theoretical framework (here iDiema nd ct al. 2005, is used) to 
interpret the spatial bias of the MPGCs, relative to the host 
galaxy's mass profile. MPGC kinematic information from 
line-of-sight velocity measurements are an optional 4th in- 
gredient that can be used for an additional, albeit weaker, 
constraint on ua. The kinematics also prove to be a useful 
tool to constrain the galaxy mass model in a self-consistent 
manner to the theoretical framework. 

The main physical model that the following analysis 
depends upon is outlined here: 

• Each MPGC formed at high redshift in a small, rare 
dark matter halo (likely together with other MPGCs and a 
small proto-galax;y). 

• The bath of UV radiation associated with the reioniza- 
tion epoch was such that MPGC formation was "instanta- 
neously" suppressec{3 in the GC formation sites. This allows 
one make a simple link between the MPGCs and a single i/a 
class of halos, thereby constraining Zreion exactly. 

• The halo containing MPGCs was later accreted by the 
host galaxy halo, i.e. it became a subhalo. 

• This accretion process effectively disrupts the early ac- 
creted, biased, rare dark matter subhalos, but the globular 
cluster tracers remain intact and will maintain the spatial 
and kinematic properties of the subhalos even to the present 
day (Diemand ct al. 2005). 

Should the second assumption prove to be incorrect, 
the analysis still provides an interesting characterisation of 
the progenitor subhalos of MPGCs and perhaps the other 
mechanism(s) that truncated their formation. Also, the trac- 
ers need not be MPGCs and can include e.g. the Milky Way 
stellar halo (IDiemand et al. 2005.) and its satellite galaxies 
iMoore et al.ll2006l. see also lOcvirk fc Aubertll201l[ ). Indeed, 
a useful feature of the Diemand-Moore technique is that it 
does not care about the detailed baryon physics of e.g. star 



^ If a locally-instantaneous reionization model is incorrect and 
an extended reionization is preferred, the MPGC constraints still 
provide a lower-limit on zreion- This is because the MPGCs will 
tend to be associated with the most abundant class of halo at 
the end of the extended reionization epoch: the least-massive or 
lowest fa halo that has only just exceeded Tvir = 10* K and 
can form MPGCs. The final spatial distribution of MPGCs will 
therefore reflect this low ua class of halos and hence the epoch 
when reionization completed can still be recovered. 
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formation, but only requires an observable tracer of subhalos 
that accreted relatively early in a galaxy's formation history. 



3 TARGET GALAXIES AND THE 
OBSERVATIONS 

In this Section, the observational data are described. To con- 
strain va for MPGC progenitors in a galaxy, the MPGC 
surface density profile must be measured out to large radii 
(typically tens of arcminutes), where contamination starts 
to contribute significantly to the profile. To reduce contami- 
nation levels, high-quality imaging is needed. To statistically 
subtract the remaining c ontamination, spect roscopic infor- 
mation can be used (e.g. IStrader et al.|[201]}) or the profile 
must extend to large enough radii to characterise contami- 
nation levels. Furthermore, the galaxy must have relatively 
"well-behaved" kinematics so that an estimate of the mass 
model can be derived. 

The extragalactic targets described here are the only 
galaxies that currently meet these demanding requirements. 
The galaxies also conveniently reside within different galaxy 
environments, including the "high-density" environment of 
a cluster galaxy, a "medium-density" group and a "low- 
density" field galaxy. 



3.1 Messier 87 

Messier 87 (M87) is the central galaxy associated with the 
Virgo cluster of galaxies. It is located 16.5 Mpc away and the 
cluster contains hundreds of member galaxies. Though M87 
and the Virgo cluster are frequently treated as the same 
thing (at least in terms of their "common" dark matter 
halo), recent analysis supports early suggestions that the 
Virgo cluster is not a single massive galaxy cluster halo, but 
is instead made up o f a set of subhalos in an unrelaxed state 
l|Binggeh et al.ll 19871: [Strader et al.ll201ll) . The M87 galaxy- 
sized halo is thought to be virialized, but is a distinct sub- 
halo from the rest of the Virgo structure (see discussion in 
Dohertv et al. I l2009l : IStrader et allboill: iRomanowskv et all 



20 111) . 



From kinematic tracers within ~ 100 kpc, the ex- 
trapolate d virial mass (Myjr ) of the M87 halo is Afvir 



10" Mq (|Strader et allboill ). The surrounding Vi rgo clus- 
ter h alo is Mvi. ^ lO^'^ Mrr^fe.g. lMcLaughlinlll999bl: Strader 
et al.lioi^Twhicii means M87 resides in high-density en- 
vironment and was therefore one of the first locations in 
the local Universe to host star formation. It has even been 
proposed that ionizing radiation from sources in Virgo had 
enough of a head-start to travel t o the Local Group an d 
reionize the Milky Way galaxy (e.g. IWeinmann et al]|2007h . 
For these reasons, M87 and the Virgo cluster are a partic- 
ularly interesting location to measure the local reionization 
epoch. 

The M87 GC catalogue is presented in IStrader et al.l 

(2011) and was constructed from archival CFHT/Megacam 
(Bouladc ct al. 20 03) images. Kinematic information for 



the GCs come from IStrader et al.l (|201ll ). who p resent new 



Keck /DEIMOS (|Faber et al.ll2003l), Keck/LR IS (jOke et al.! 
[liii) and MMT/Hectospec l|Fabricant et al.| [2005') hnc-of- 
sight velocities. Out of an estimated MPGC population of 



~ 14000 MPGCs l|Strader et al.ll201ll ). 289 have useful ve- 
locity information. 

The M87 MPGC surface densitv was derived in Strader 
et al. pOllf) using standard reduction techniques and 
methodologies. Briefiy, GCs were photometrically selected 
down to i = 22.5 mag. over a region extending to radii 
~ 130 kpc from M87. Hubble Space Telescope imaging and 
the velocity information were used to accurately constrain 
the amount of contamination present in the GC catalogue, 
which was subtracted off the MPGC surface density profile 
before analysis. 



3.2 NGC 1407 

The second target of this study is an elliptical galaxy named 
NGC 1407, whose g alaxy group contains ~ 30 members 
(|Brough et al.l |2006|). It is 21 Mpc d i stant and has a siz - 
able GC system llForbes et al.l [200^ : iHarris et ai1l2006l ). 
IRomanowskv et all ("200^ found the galaxy group, which is 
centred on NGC 1407, to have a relatively high mass-to- light 
ratio for its total virial mass Afvir ~ 6 x 10^"^ Mq. This tar- 
get provides a probe of reionization in a "medium-density" 

environment. 

Two Subaru/Suprime-Cam (|Mivazaki et al.ll2002l ) im- 
age datasets of NGC 1407 were analysed. Each covers a 
~ 34' X 34' area on the sky. The first is centered on 
NGC 1407 and is made up nearly 5 hours of g, r, i imag- 
ing under 0.5 — 0.6" conditions. The resulting GC cata- 
logue has been used fo r an extensive spectroscopic campaign 
with Keck/DEIMOS (IRomanowskv et al.ll2009l : iFoster et all 
I2OIOI: IPota et aLlbOllT and a photometric studv ( Forbes 
et al. l201l[) . The second Suprime-Cam dataset is an archival 
dataset in B, V, I bands with 0.6 — 0.7" seeing and totaling 
~ 2 hours in exposure time. In this imageset, NGC 1407 was 
positioned toward one end of the mosaic so the radial cov- 
erage extends ~ 30' from NGC 1407 or ~ 200 kpc. The GC 
catalogues are available upon request from the first author. 

MPGC surface density profiles were constructed in each 
dataset from a catalogue of GC candidates found in the im- 
ages (MPGCs were taken to have (g — i)o ^ 0.98). A liter- 
ature Hubble Space Telescope, A dvanced Camera for Sur- 
veys catalogue (|Forbes et al.ll20()^ was incorporated into the 
analysis to improve the surface density profile in the central 
~ 3' of the galaxy center. Standard procedures were fol- 
lowe d to derive the GC surface densitv profile (e.g. Spitler 
et al. |2008bf). The surface density data extend to ~ 210 kpc 
and a constant surface density is not reached at large radii, 
suggesting a small number of MPGCs should be found be- 
yond the region covered by the Suprime-Cam imaging. Since 
contamination was not subtracted from the MPGC surface 
density profile, the fits to the NGC 1407 MPGCs described 
in subsequent sections include a constant background level 
to model the contamination. 

An expanded GC lin e-of-s ight velocity catalo gue from 
IRomanowskv et al.l (|2009l) and iFoster et al] (|2010l) is used, 
whic h incorporates new Keck/DEIMOS observations ( Pot a 



et al. I2OIII) . The total number of NGC 1407 MPGCs with 
spectroscopic velocities is 167 GCs (excluding suspected ul- 
tra compact dwarf candidates; see details below), out of an 
estimated total of ~ 4000. 
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3.3 Milky Way 

The Milky Way (MW) MPG C spat i al and kinematic infor- 
mation come from the new iHarrid \20ldb catalogue. The 
MPGC sample contains 110 GCs with [Fe/H] < -0.9 and 
has a mean metallicity of [Fe/H] ~ —1.6. The Milky Way 
GC analysis is carried out in three-dimensions rather than in 
projection as for the preceding gala^des. In order to model 
the halo mass profile of the Milky Way with the MPGC 
velocity information, heliocentric distances and velocities 
must be corrected to the Galactic reference frame. The dis- 
ta nce correction is straightforward and included in the Har- 
ris l|2010ll catalog itself, but the velocities in principle re- 
quire proper motion measurements, and these are in many 
cases not available with sufficient precision. Instead, after 
correcting the line- of-sight velocitie s for the heliocentric mo- 
tion (e.g., eq. 5 of IXue et al.ll200^ . the data are corrected 
to the true radial velocities on a statistical basis: given the 
anisotropy profile P{r) from the best-fit model (discussed in 
the next Section), the corre ction from Une-of-sig ht to radial 
velocity dispersion (eq. 1 of lBattaglia et aL l liooQ) is applied 
to the observed velocities, and the dispersion profile is calcu- 
lated. This processes is repeated iteratively and ultimately 
yields a very small correction that affects the final u values 
at only the ~ 3% level (after excluding from the analysis the 
GC kinematics inside r = 8 kpc, which are more difficult to 
correct). 



4 ANALYSIS 

The Diemand-Moore technique uses the spatial bias of 
MPGCs relative to the hosting galaxy's dark matter dis- 
tribution to constrain the i/a properties of the MPGC pro- 
genitor subhal o s. Th is is accomplished by fi tting a modified 
iNavarro et al.1 l| 19961 . NFW) function from iDiemand et al.1 
( 2005h to the MPGC spatial distribution. The fit also re- 
quires as input the scale radius (rs) of the host galaxy halo, 
and thus an accurate mass model of each galaxy is an im- 
portant component of this work. Although models for these 
galaxies exist in the literature, they are derived using differ- 
ent techniques and various mass tracers. To remove possi- 
ble systematic errors resulting from these inhomogeneities, 
new mass models are derived for each galaxy. Also, to en- 
sure the ua fits to the surface density profile are derived in a 
self-consistent manner with the galaxy mass models, MPGC 
observations are used to derive the mass model. 

The mass modelling methodology is outlined in the first 
subsection below. In subsequent subsections, the individual 
galaxy mass models are derived. Though some necessary as- 
sumptions and compromises are made to yield the final mass 
models, the halo mass properties and their uncertainties are 
robust. In the final subsection, ij4.6l the measurements 
for each are converted into Zreion constraints. 

In the following, the MPGC spatial profile is denoted 
j{r) and S(_R), for three-dimensional and the projected 
cases, respectively. Also, when ua constraints from the 
Diemand-Moore technique are discussed, f actually refers 
to all subhalos above the subhalo virial mass correspond- 
ing to that I'. For example, the constraint v = 3.1 for the 
Milky Way, actually includes all subhalos equal to and more 
massive than a, u — 3.1 halo at the redshift of interest. 



4.1 Derivation of the Galaxy Mass Models 

Here we outline our modeling methods that will be used for 
specific galaxies in the following subsections. We begin with 
the three-dimensional number density profile of the dynami- 
cal tracer population (MPGCs in this paper). This profile is 
mod eled as a modified NFW profile, as derived in Diemand 
et al. (|2005^ : 

"'^ = (r/r-.)-r[l + (rV)'^](''--^)/'^' 

where u is the rarity value, 7 and /3„ are the inner and outer 
power-law slopes, respective ly, and a descr i bes th e inter- 
mediate behavior. Following IPiemand et aP ()2005l ). we set 
Q = 1, 7 = 1.2, and = 3-|-0.26!/^-®. Thus a higher u value 
implies a steeper outer slope. The value of v is related to a 
characteristic radius ri, by: 

u = 2\n{rs/r^), (2) 

which means the surface density is more compact for higher 
values of u. 

Ideally, we would compare this model to an observed 
surface density profile, and fit for 1/. However, in practice 
there are not good prior constraints on Vs available, i.e., we 
do not know the DM distribution of the halo which we need 
to estimate the relative compactness of the MPGC subpop- 
ulation0 To solve this problem, we will make use of the fact 
that there are kmematical predictions of the model, so we 
can use the observed spatial and kinematical distributions 
to fit simultaneously for v and r^. 

Our basic approach for the kinematics is to take model 
distributions for mass, MPGC d ensity, and anisotropy, solve 
a spherical Jeans equation fe.g.. lMamon fc LokadbOOai and 
derive a velocity dispersion profile for the MPGCs, cr(r). In 
the Milky Way, the (j(r) prediction can be directly com- 
pared to observational data, while in external galaxies, an 
additional step is needed to project to line-of-sight disper- 
sions, ap{R). 

Each mass model consists of a stellar component (with a 
reasonable value for the stellar mass-to-light ratio T«) and a 
ACDM-based model for the DM distribution. This DM pro- 
file has a similar parametric form to Equation [1] but with 
density and scale parameters {ps and r^) that are highly cor- 
related. To derive these correlations (or modeling priors), we 
begin with the relation between virial mass and concentra- 
tion predicted for relaxed halos at z = by the Bolshoi and 
BigBo lshoi simulations u sing WMAP5 cosmological param- 
eters ()Prada et al.ll201ll , eq. 25). After converting the mass 
and concentration into scale radius and density, and sub- 
tracting off a cosmological baryon fraction of 0.17 to find 
the DM density alone, the following power-law provides an 



^ One might think that with massive ellipticals like NGC 1407 
and M87, an X-ray based mass profile would provided the much- 
needed inde pendent information on t he DM profile. However, as 
discussed in lRomanowskv et al.l ||2009|) . for these particular galax- 
ies, there are inconsistencies between the X-ray results and GC 
system dynamics. In fact, there is an emerging, wider pattern of 
disagreement between X-ray and optical mass determinations in 
elliptical galaxies (e.g.. [Shen fc Gcbhardt 2010). For this paper, 
it is assumed that the MPGC dynamics results are reliable. 
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Figure 2. Modeling ingredients for the M87 metal-poor globular cluster system. Model predictions are shown as curves: the solid curves 
are the best model, with u = 2.19, Ts = 126 kpc, and ps = 1.4 X 10^ Mq kpc~^ (a normal concentration halo with Mvir ~ 4.8 X 10^'^ Mq). 
The dotted curves are u = 2.61, = 191 kpc, and ps = 0.7 X 10® Mq kpc"^ (a low concentration halo with Mvir ~ 7.2 X W^^Mq); and 
u = 1.80, Ts = 85 kpc, and pa = 2.9 X 10® Mq kpc~^ (a high concentration halo with Mvir ~ 3.4 X IO^^Mq). Note, the mass modelling 
only incorporates tra cers from within ^ 100 kpc of M87 so th e halo properties are a reflection of the M87 subhalo, not the entire Virgo 
cluster halo (e.g. see lDohertv et aljlioosl: IStrader et alj|201ll) . R is for projected radii, r is for physical 3D radii. Upper left: Circular 
velocity profile. The total profile for the best model is shown decompose d into its stellar and dark matter co mponents (dashed curves, 
as labeled). Independent results are also shown from stellar dynamics jMurphv et alj|201lf) and X-ray gas llDas et alj|201oh . Upper 
right: MPGC velocity anisotropy in the model (lines) and estimated from the kurtosis of the observed line-of-sight velocity distribution 
(circl es). Lower left: MP GC surface density. Lower right: MPGC line-of-sight velocity dispersion. Data are shown by points with error 
bars jStrader et al.|[201lh . In fitting the dispersion data, there is a strong degeneracy between the halo parameters and ps, which is 
broken only by imposing cosmological priors. The remaining uncertainty on translates directly into an uncertainty on u. The dashed 
curves show two models that are ruled out by the dispersion data, even after allowing for a 1 cr scatter in the concentration. 



approximation to the average results of the numerical sim- 
ulations; 

which is good to 8% or better over the virial mass range 
of ~ 10^^-10^* Mq. a 1 a cosmological scatter in halo con- 
centr ation of ±0.11 dex (iDuffv et al.ll2008l: iMaccio et al] 
l2008h translates to a scatter in ps at fixed Vs of ±0.24 dex. 
This scatter ultimately provides the largest systematic un- 
certainty on the measurements of i/ and hence Zreion- 



These results involve fits of standard NFW density pro- 
files to the simulatio ns (i.e., Equati on [1] with 7 = 1, a = 1, 
P = 3). However, the lPiemand et al.l ()2005l) analyses used a 
slightly different DM density model, where the central cusp 
has a power-law slope of 7 = —1.2 rather than —1: 

Pir)^p.[^y\i+^y\ (4) 

Here it is assumed if such a model had been fitted to the 
recent simulations, the Afvir and values would be the same 
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as with the NFW fits. Eq. [3]then becomes: 



1.07 X 10' 



kpc-^ 



kpc 



(5) 



In practice, the increased fractional mass inside means 
that lower ps values will be fitted to the data, relative to an 
NFW model. The model for a mass-sequence of average DM 
circular velocity profiles is then: 



i'c,DM(r) = 17.9 



km / rs V ™ f r 



s \kpc J 

2F1 f 1.8, 1.8; 2.8; 



11/2 



(6) 



where 2-F1 is a hypergeometric function. 

The final ingredient in the models is the velocity disper- 
sion anisotropy profile, which is normally a major source of 
uncertainty in modeling observations of pressure-supported 
galaxies, but in the context of the Diemand-Moore models, 
is uniquely constrained. This anisotropy profile is: 



where: 



/3o = 1- 



r + ri,' 



5^' 



(7) 



(8) 



This anisotropy profile is isotropic in the center (/3 ~ for 
r <^ r^) and becomes radially biased in the outer regions 
(/? ~ Po for r 3> r,^). Higher values of 1/ imply stronger 
radial anisotropy. 

A complication for modeling both NGC 1407 and M87 
is that there are some correlations between GC luminosities 
and kinematics that are probably driven by populations of 
ultra-compact dwarfs (UCDs) at the bright end, with dis- 
tinct dynamics fr om the "n ormal" GCs (see Romanowsky 
et al. 2009; Strad er et al.l l2011: Brodie et al. 2011, for more 
details). Although the UCDs comprise only a tiny part of 
the overall "GC" systems, they represent a substantial frac- 
tion of the luminosity-based spectroscopic samples. Such ob- 
jects were therefore eliminated from the kinematic samples 
by excluding anything that has a large measured size (half- 
light radius > 5 pc), or that has no size measurement but 
is bright (io < 21.3 and io < 20 for NGC 1407 and M87, re- 
spectively). These objects have no impact on the measured 
surface density profiles. 

The following subsections present the details of the 
modeling for each galaxy, and the resulting 1/ constraints. 



4.2 Rarity of M87 MPGC progenitor subhalos 

For M87, the mass models are tuned not only to fit the 
MPGC subpopulation velocity dispersion but also to have 
a total Vc — 500 km s~^ at r ~ 5 kpc for cons i stency with 
the stellar dynamics results of iMurphv et all ()201lh : this 
mea ns adiusting T„ accordinglv. As discussed in Strader 
et al. (2011), the NFW models have some difficulty in re- 
produc ing the stellar dynamic s inferences at r ~ 15 kpc; the 
cuspier iDiemand et~all (|2005l ) models fare better, but there 
may still be some tension with the stellar data. These mass 
models are also comp letely unab l e to r eproduce the X-ray 
based mass- model of iDas et all l|2010l) while maintaining 



strongly radial anisotropy and reproducing the observed ve- 
locity dispersions (an additional example of the problem is 
footnoted in H4.1|) . 

Figure [2] shows the elements of the modeling procedure, 
along with comparisons to data as applicable. The model 
consists of profiles of Vc{r), P{r), and j{r) (the latter is not 
shown), which are then input to the Jeans equation to solve 
for (7,,(r) and cre(r), and projected to E(i?) (lower left panel 
in Fig. [5| and Vrma{R) (lower right panel). As a consistency 
check, the model expectations for /3(r) are compared with 
an empirical estimate from the line-of-sight velocity distri- 
bution of the data in the upper right panel of Fig. [21 This 
estimate makes use of a restricted class of spherical constant- 
P models which predict the line-of-sight velocity kurtosis 
if the veloci ty dispersion profile is approximately constant 
with radius (|Napolitano et al]l2009^ . The kurtosis measure- 
ments from the MPGC data can thus be converted to rough 
estimates of the anisotropy profile. 

For M87, a typical-concentration ACDM halo with rs = 
126 kpc (M yjr ^ 5 X 10"Mf^) is found, along with the Die- 
mand et al.^2005^ anisotropy prediction, to fit the M87 
MPGC dispersion data remarkably well. The implied /3(r) 
is also nicely consistent with the empirical estimates. This 
model corresponds io u = 2.19. The predicted velocity dis- 
persions are fairly sensitive to the value of , increasing for 
larger rs (and v). Using a fit to the dispersion data, is 
determined at the ±10 kpc level, and u at the ±0.08 level. 

These results hold for the strong theoretical ACDM 
prior adopted on the halo parameters. Relax;ing this prior 
allows for a much larger range of solutions, as various com- 
binations of halo density and scale radius can all provide 
similar mass profiles in the region with MPGC observations. 
Allowing for the predicted 1 a scatter in the halo properties 
(i.e., varying the normalization of Eq.[6]by ± 0.12 dex), rs is 
constrained at the ~ 50% level, yielding v = 2.19±0.41 (see 
Figure [2}. 



4.3 Rarity of NGC 1407 MPGC progenitor 
subhalos 

For NGC 1407, several priors are invoked: the mass- 
concentration relation from ACDM; an estimate of the mass 
at i? ~ 100 arcmin from satellite galaxy dyn amics {vc ~ 
700 ±100 km s~^: see lRomanowskv et al.l2009t) : and a mea- 
surement of the central s tellar velocity disp ersion around 
1 Rc of cTp ~ 260 km s"^ (|Proctor et al.ll2009t ). It is beyond 
the scope of this paper to carry out the detailed dynami- 
cal modeling required to convert this ap measurement into 
an estimate of Vc , and instead the empirical results from dy- 
namical modeling of vario us other pressure-s upported galax- 
ies t hat vJar. ~ 1.4-1.7 dWolf et al.ll2010l: Truii llo-Gomez 
et al. l201ll : lMurp° et al .IbOlltlDutton et al ] |201lh are used. 

Unlike M87, for NGC 1407 there is no model solution 
that straightforwardly fits all of the above constraints (see 
Figure [3]). In particular, it appears to be very difficult to 
match the inner MPGC velocity dispersion profile, where 
the CTp ~ 200 km s"'^ inside ~ 4 arcmin (~ 25 kpc) is lower 
than the stellar velocity dispersion. Normally, the MPGC 
dispersion should be higher because the spatial distribution 
is more extended. Lowering the dispersion would require tan- 
gen tial orbital anisotropy , which contradicts expectations in 
the lDiemand et al] (|2005l ) model for mild radial anisotropy. 
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It is speculated instead that GC disruption processes near 
the galactic cent er have preferentially depleted the more ra - 
dial orbits (e.g. iFall fc Zhan3l200ll : FVesperini et al.ll2003h . 
MPGCs inside 3.5 arcmin (~ 20 kpc) are therefore excluded 
from the analysis (it is not clear why this should be neces- 
sary for NGC 1407 but not for M87). 

Considering now the outer MPGC velocity dispersion 
data, it also turns out to be difficult to fit the joint con- 
straints from ACDM halo concentration expectations and 
the satellite galaxies dynamics. The mass required by the 
latter predicts MPGC dispersions higher than observed, 
which can be alleviated by a low-concentration halo. The 
range of plausible solutions is bracketed by adopting a halo 
that has a low concentration at the 1 o level while fitting 
the satellites constraint, or one with a normal concentration 
that does not agree with the satellites (whose constraint is 
then assumed to be unreliable). Alternative modifications to 
lower the MPGC dispersions would be to use a less cuspy 
NFW model (which is within the cosmological scatter for 
halo profiles), or to decrease the adopted distance since the 
kinematics data will then effectively be probing less far into 
the halo (although this galaxy is likely to be more distant, 
not less, based on estimates summarized in NED). 

The resulting range in fit parameters is z/ ~ 2.6-3.1; 
again, the systematic uncertainties outweigh the statistical 
ones (~ ± 0.15) from fitting the dispersion and density pro- 
files. An additional issue from these solutions is that the low 
observed kurtosis for the outer MPGCs suggests tangential 
anisotropy in these regions too, again at odds with the radial 
model expectations. 



4.4 Rarity of Milky Way MPGC progenitor 
subhalos 

The Milky Way modeling differs slightly from the previous 
two galaxies by treating 3-D rather than projected profiles 
(of density and velocity dispersion). For the mass model, a 
ACDM motivated halo is included as before, along with stel- 
lar d isk and bulge compon ents based on models "I" and "H" 
from lBinnev fc Tremain^ ((2008i) . which bracket a plausible 
range of parameters for th ese component s . We also tried 
the disk and bulge model of IWidrow et al.l (|2008l l. which is 
roughly halfway between the I and II models of Binnev & 
Tremaine (|2008h . As shown in the upper left panel of Fig. U 
since these components influence the central regions only, all 
choices give the same results for v. The model I of Binnev 
& Tremaine (|2008l ) was adopted. 

The modeling results are illustrated in Figure U After 
excluding the region inside r = 2 kpc f rom the density proflle 
fit because of the pres ence of a core (|Parmentier fc Grebell 
l2005l: iBica et al.ll200d) . a good fit to both the density and 
radial velocity dispersion profiles is found. The best-fit value 
for the MPGC subhalo rarity is = 3.12 ± 0.34, where the 
uncertainty is again driven by the scatter in the assumed 
halo mass-concentration relation. The parameters for the 
DM halos in these models are listed in the Figure caption. 
The different DM model derived here is the main reason why 
a higher i/ is foun _ d compared to the earlier work bv Moore 
et al. l200tl . 

Note that there is a gap in the MPGC velocity data 
around r ~ 50 kpc. If the MPGC velocity dispersion in this 
region were assumed to be as low as for the halo stars (as 




10 15 
redshift z 



Figure 5. Evolution of subhalo rarity {va) over a range of virial 
masses and redshifts z. Various va peaks are indicated by dotted, 
grey lines. Solid black line shows the mass-redshift pairs, where 
halos have virial temperatures Tvir = 10"^ K and efficien t star 
formation can take place (curve from iMoore et al.ll2006ft . The 
constraints for the galaxies studied here are given as the coloured, 
solid lines: v = 2.2± 0.4, v = 2.9 ± 0.4, and u = 3.1± 0.3 from 
left to right (M87, NGC 1407 and the Milky Way, respectively). 



indic ated bv the dashed curve in the Figure, from Gnedin 



et al. |201G| ) then the virial mass and v value would need to 
be lower. The default ha l o mas s profile is actually similar 
to that of iGnedin et al.l l|2010l) despite the higher veloc- 
ity dispersion found here. This is because the high typical 
anisotropy j3 ~ 0.65 in the current model elevates the radial 
dispersion relative to the (3 ~ 0.4 that the former authors 
assumed, and which would be appropriate for the dark mat- 
ter particles on average, not for the "biased" subset of halo 
particles that should correspond to the stars and MPGCs. 

4.5 Mass modelling summary 

In the Milky Way and M87. the MPGC k inematics comfort- 
ably fit within the lDiemand et aL 1 (|2005D theoretical frame- 
work and are consistent with other mass tracers. Though 
the mass models unfortunately suffer from rather large sys- 
tematic uncertainties from the ACDM prior, the quality of 
the models mean the va measurements described below are 
on firm ground. The situation for NGC 1407, from a mass- 
modelling perspective, is less clear. There appears to be lin- 
gering tension between the models and observations, which 
merit further study. 

4.6 Constraining the reionization epoch 

To convert the uo constraints from the preceding sec- 
tions into estimates for the truncation redshift of MPGC 
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Figure 3. As Figure [2] for NGC 1407. Error bars in tiie upper-left panel are mass profile constraints from stellar and satellite galaxy 
dynamics, as discussed in the text. Two possible solutions are shown by solid curves: one is a normal-concentration halo with i/ = 2.61, 
= 122 kpc, ps = 1.4 X 10^ Mq kpc-3, and M^ir ~ 4.4 X IO^^Mq; the other is a low-concentration halo with u = 3.11, = 205 kpc, 
Ps = 0.7 X 10'^ Mq kpc-3, and M^ir ~ 8.6 X 10^^ Mq. 



Table 1. Rarity (ua) of MPGC progenitor subhalos, the reion- 
ization redshifts and minimum virial mass of the progenitors for 
the target galaxies. 





environmental 
density 


V 




min. virial 
mass 
[/i-i Mq] 


Messier 87 
NGC 1407 
Milky Way 


high/cluster 
medium / group 
low/field 


2.19 ±0.41 
2.86 ±0.36 
3.12 ±0.34 


7.8" 
11.0' 
12.1" 


f2.0 
-1.8 
fl.7 
-1.7 
H.6 
-1.1 


1 X 10* 
8 X 10^ 
7 X lO'' 


mean 




2.74 ±0.21 


10.5" 


n.o 

-0.9 


8 X lO'^ 



formation (e.g. Zreion), a degeneracy must be broken be- 
tween subhalo virial mass and redshift for a given v. The 
evolving properties of a gi ven va subhalo can be analyti- 
cally computed using the (|Press fc Schechteilll974l ') Gaus- 
sian field formalism. The value v is explicitly defined as: 



V = (5c/[i?(2:)(7(Mvir)], where the constant 5c is the crite- 
rion for spherical collapse, D{z) is the linear growth factor 
and (7(Mvir) is the typical mass-density fluctuation "height" 
for a subhalo of a given virial mass M^,-„ fsee e.g. Mo & 
White l2002i r For the present work, a top hat filter is used 
to select the characteristic spatial scale. The CAMB ( Lewis 
et al. JOOO) online intcrfactjf] was used to ge nerate the linear 
powe r spectrum for a WMAP7 cosmology jKomatsu et al.l 
boilh . Figure [5] shows the relationship between the rarity, v, 
of a mass-density fluctuation for a given Mvir— redshift pair. 

Figure [S] also shows how a vg measurement alone will 
not constrain the redshift when MPGC formation is sup- 
pressed. Further information i s requ ired to break this degen- 
eracy. Following iMoore et al. l l|2006l) . it is assumed that star 
and MPGC formation can only occur in subhalos that have 
virial temperatures greater than Tvir = 10*K, where gas is 



http : //lambda.gsfc .nasa.gov/toolbox/tb_camb_f orm. cfm 
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Figure 4. As Figure [J] for the Milky Way. Here tiie lower panels are three-dimensional rather than p ro jected quant i ties. T he dashed 
curve in the lower-right panel shows the velocity dispersion profile fit to halo tracers (mostly stars) from iGnedin et al.l 1I2OIOI ) . The solid 
curve in this panel shows the best model, with v = 3.12, = 33 kpc, and ps = 2.5 X Vfi Mq kpc"'' (a normal concentration halo with 
M^ir ~ 1.8 X IO^^Mq).. The dotted curves show a high-concentration halo with v = 2.79, Ts = 22 kpc, ps = 5.1 X 10^ Mq kpc ^, and 
Mvir ~ 1.3 X IO^^Mq; and a low-concentration halo with v = 3.48, = 49 kpc, ps = 1.2 X IO'^Mq kpc'^, and Mvir ~ 2.5 X lO^^M©. 



able to rapidly cool and form stars (e.g. lOstriker fc Gnedinl 
Il996h . The region in Fig.[5]below the Tvir = curve can 

therefore be excluded. The region above the T^ir = ICK 
curve can be excluded as well, since MPGC formation would 
have taken place in subhalos that would ultimately have less 
or more concentrated spatial distributions than what is ob- 
served in the galaxies. For example, for the Milky Way con- 
straint of !/ = 3.1, at redshift of 2 = 10, subhalos with mass 
Mvir ~ W^Mq could still have formed MPGCs efficiently, 
which would have led to a spatial distribution resembling 
2.5a subhalos, contradicting the observations. At redshifts 
higher than the intersection between the Tvir = 10*K rela- 
tion and the curve from MPGC observations, only rarer, 
e.g. u > 3.1, subhalos could form MPGCs, thus the spa- 
tial distribution would have been more concentrated than 
is observed. The intersection between the ua constraint and 
the Tvir = 10*K curve is therefore taken to be the MPGC 
truncation redshift or the local Zreion of the host galaxy. 



The MPGC ua fits and the corresponding reionization 
redshifts are given in Table [T] for each galaxy. The weighted 
mean of the three vmpgc measurements translates into a 
reionization redshift of Zraion = lO.SlJ g. These constraints 
are discussed in Section [S] below. 



5 DISCUSSION 

The reionization redshift (zreion) constraints from metal- 
poor globular cluster observations (see Fig. [6]) are summa- 
rized in Fig. [7] with existing estimates and limits from the 
literature. The local constraints from metal-poor globular 
cluster (MPGC) observations are tabulated in Tab. [1] 

The only existing constraints on Zreion from MPGCs 
are for the Milky Way. As described in TO our estimate 
agrees with that derived in I Moore et al.l 1^06), who used 
the same technique employed here. Also, the high-resolution 
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Figure 7. Reionization constraints from the current work and the literature. The reionization epoch measurements for the local metal- 
poor globular cluster observations (coloured constraints) are labeled with their local environments (the Milky Way, NGC 1407 and M87 
for low-, medium- and high-density environments, respectively) and dominated by un certainties on galaxy's mass model. High-rcdshift 
constrai nts include: co smic microwave background result from 'Komatsu ct al.' (2011") , constraints from UV galaxy observa tions from 
[Roberts on et al.l | |2010|) . Lyman-a emitters from Kashikawa et al. (2011); Pcntcricci ct al. (2011), and Lyman-a forest from lF^n et all 
Arrows indicate an estimated limit on the reionization epoch, while each hatched region contains the 68% uncertainty about a 
measurement. 



dark m atter Milky Way simulation described in lGriffen et al.l 
(|2010l) can only reproduce the observed spatial distribution 
of MPGCs and their numbers only when MPGC formation 
was suppressed by reionization at Zreion ~ 13. This value 
agrees with our own estimate within its la uncertainty. 

5.1 Environmental propagation of reionization 

The MPGC Zreion values and their uncertainties span the 
redshift range covered by Zreion estimates from the liter- 
ature, with the weighted mean MPGC constraint show- 
ing best agreement with the cosmic microwave background 
(CMB) v alue and an instantaneous reionization model ( Lar- 
son et al. I2OIII : iKomatsu et al.ll2011^ . It is apparent that 
the M87 Zreion measurement is somewhat of an outlier com- 
pared to the other MPGC estimates. Since M87 is located in 
the high-density environment of a galaxy cluster, it is pos- 
sible that its lower Zreion results from an environmentally- 
dependent reionization epoch. Unfortunately, the relatively 
large systematic uncertainties on the galaxy halo mass mod- 
els mean the difference between Zreion in the highest and 



lowest density environments (M87 vs. Milky Way Zreion) is 
only significant at the 1.7a level. 

In addition to a correlation between the MPGC spa- 
tial distribution and the local reionization epoch. Moore 
et al. (|2006l ) predicted that Zr-^inr, should also correlate 
with the virial mass-normalized MPGC numbers. This is be- 
cause gala^xies with a lower Z r eion ha d mor e time to produce 
MPGCs. Both ISpitler et all (|2008bl ) and ISpitler fc Forbes! 
( 2 009. ) found no evidence for a correlation between MPGC 
virial mass specific frequencies and environment. However, 
their virial mass estimates were derived from statistical re- 
lationships between stellar and virial masses, which are not 
ideal for use on individual galaxies. 

Using the galaxy virial masses derived in Section [4] and 
MPGC GC numbers (Nmpgc) from Sectional the follow- 
ing MPGC virial mass specif ic frequencies (VMpan Spitler 
et al. 2008h) are found: Vmpgc = Nmpgc /{M^ir/W^'^) = 
27.5l8^/,6.7li;9,6.lli;7, for M87, NGC 1407, and the 
Milky Way, respectively. The errors are again dominated 
by systematic uncertainties in the galaxy mass models. M87 
has a Vmpgc value that is larger than the Milky Way or 
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Figure 6. Normalized, projected surface density profiles of the 
metal-poor globular clusters in each of the target galaxies, which 
are labelled in the figure. Curves show the subhalo spatial dis- 
tribution for a range of rarity values (e.g. v ^ 2, 3, 4 standard 
deviations from the cosmic mean mass-density), which directly 
constrain the epoch at which metal-poor globula r cluster forma- 
tion was halted by the local reionization epoch llDiemand et alj 
I2OO5I : iMoore et alj200^ .From fitting the theoretical expectations 
for the spatial and kinematic distributions of metal-poor globular 
clusters (see the epoch of reionization is constrained to be 



12.1^ 



-,-1-1.7 



:^for 



the Milky Way, NGC 1407 and M87, respectively 



NGC 1407, at the 2. So" significance level. The normalized 
MPGC numbers thus lend further support to the idea that 
reionization finished in M87 at lower redshifts compared to 
the other galaxies. In general, massive galaxies at the centres 
of galaxy clusters (e.g. M87) are known to have enhanced 
relative number of GCs when normalized to th e galaxy's 
star light (e.g. lHarrislll99ll: lMcLaughlinlll999ah . Whether 
this holds when the MPGC numbers are instead normalized 
to virial masses (which should not be directly influenced by 
reionization) will need to be explored in a future work. 

The above results provide an exciting hint that reion- 
ization completed first in low-density environments. 



5.2 The theoretical context 

A number of large-scale radiative transfer simulations have 
been run to understand the propagation of reionization 
through different galaxy environments. From analysis of 
these simulations, it is found that proto-galaxies in high- 
density environments produce significant numbers of ioniz- 
ing sources early on and thus are the first l ocations to reion- 
ize (llliev et al.ll2006l: iFinlator et al.|[2009l: llliev et~aLll201ll 



dynamical range to capture both the environmental prop- 
agation of reionization over large volumes and the small- 
scale physics that are important for understanding localized 
affects. For example, when an inhomogeneous intergalac- 
tic medium is incorporated into reionization models, dense 
pockets of neutral gas can survive even after surrounding 
regions had completed the reionization process (Miralda- 



Escude et al. 2000': 'Ciardi ct al."2003":'Furlanetto fc Ohll2005 



Iliev et al.i,2008. : Choudhury et al. 2009 ; Mitra et al.l 12011 
2012h . These pockets tend to be located in dense environ- 



ments, where recombination rates can stay high despite their 
proximity to large reservoirs of ionizing sources. This means 
the intergalactic medium in dense environments started to 
reionize first, but only completely reionized a/£er lower den- 
sity environments had finished. 

The trend hinted at by the metal-poor globular cluster 
observations agrees with this latter theoretical work. Under 
this scenario, MPGCs in a high-density environment contin- 
ued forming to lower redshifts in self-shielded pockets, even 
as low environmental-density regions were being completely 
ionized. The progenitor molecular clouds of MPGCs located 
in more vulnerable environments (like that surrounding the 
Milky Way) were not surrounded by gas with sufficiently 
high recombination rates to survive the ionizing front that 
originated within or passed through the region. 

Finally, the range of the MPGC constraints on Zrt^ion 
agrees with the range of literature constraints (see Fig. [7|. 
This seems to support models of reionization tha t are g lob- 
ally extended in duration (e.g. iMitra et alll201ll . 120121 ). 



however see lWeinmann et al.ll2007l ) 



Such simulations typically do not have a large enough 



6 CONCLUSIONS 

The spatial and kinematic properties of metal-poor globu- 
lar clusters (MPGCs) provide an exciting hint that reion- 
ization was inhomogeneous. In the local Universe, reion- 
ization completed first in low-density environments around 
Zreion ~ 11 — 12 (for the Milky Way and large central group 
elliptical NGC 1407) and finished last in the high-density en- 
vironment around the Virgo cluster at Zreion ~ 8. While un- 
certainties on the galaxy halo mass models limit the strength 
of this conclusion, another property, the relative number of 
metal-poor globular clusters, supports this interpretation. 
The apparent environmental dependence of Zreion also aligns 
with theoretical work that finds that the ionization front 
propagated from high- to low-density environments, with the 
high densities finishing reionization last. Furthermore, the 
apparent agreement between the range of Zre.ion constraints 
from MPGC observations and those from the literature may 
provide additional support to models where reionization was 
globally prolonged in duration. 

More observational work, especially to improve the in- 
dividual galaxy halo mass models, is required to confirm 
the results found here. It will also be important to model 
MPGC formation in large-scale, radiative transfer simula- 
tions with realistic models of star cluster dynamical evo- 
lution. This will help refine the technique utilized here and 
clarify the relationship between reionization and the trunca- 
tion of MPGC formation. For example, given that high star 
formation rates and gas densities were needed to produce 
such dense star clusters, it is possible that MPGCs trace 
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the end of reionization, where ionizing radiation was able to 
finally penetrate the densest gas clouds. 

Expanding this initial study to other galaxies will yield 
a reionization map of the local Universe, further aiding our 
understanding of galaxy formation and evolution over cos- 
mic time. 

If reionization is not responsible for the truncation of 
MPGC formation, the results presented here provide unique 
constraints on the formation histories of MPGCs. Further- 
more, since MPGCs and stellar h alos of the M i lky Way share 
spati al and chemical properties llHelmill2008l : iMartell et al.l 
I2OIII . see also M31, iHuxor et al.ll201ll) they might share a 
common origin. If this is the case for other galaxies, then 
MPGC observations can be used to constrain the star for- 
mation history of their stellar halos. 

Finally, the power of the Diemand-Moore technique is 
that it does not care about the detailed star formation and 
enrichment physics - it only requires an observable tracer 
object that accreted onto the galaxy at earlv times ( Die- 
mand et al. l2005h . 
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